Evolution of the Cluster Correlation Function

We study the evolution of the cluster correlation function and its richness-dependence from z = 0 to z = 3 using large-scale cosmological simulations. A standard flat LCDM model with \Omega_m = 0.3 and, for comparison, a tilted \Omega_m = 1 model, TSCDM, are used. The evolutionary predictions are presented in a format suitable for direct comparisons with observations. We find that the cluster correlation strength increases with redshift: high redshift clusters are clustered more strongly (in comoving scale) than low redshift clusters of the same mass. The increased correlations with redshift, in spite of the decreasing mass correlation strength, is caused by the strong increase in cluster bias with redshift: clusters represent higher density peaks of the mass distribution as the redshift increases. The richness-dependent cluster correlation function, presented as the correlation-scale versus cluster mean separation relation, R_0 - d, is found to be, remarkably, independent of redshift to z <~ 2 for LCDM and z <~ 1 for TCDM (for a fixed correlation function slope and cluster mass within a fixed comoving radius). The non-evolving R_0 - d relation implies that both the comoving clustering scale and the cluster mean separation increase with redshift for the same mass clusters so that the R_0 - d relation remains essentially unchanged. The evolution of the R_0 - d relation from z ~ 0 to z ~ 3 provides an important new tool in cosmology; it can be used to break degeneracies that exist at z ~ 0 and provide precise determination of cosmological parameters.Comment: AASTeX, 15 pages, including 5 figures, accepted version for publication in ApJ, vol.603, March 200

Dynamic Interaction Between Rail Track Systems and the Layered Subsoil Solutions in the Frequency-and Time Domain

For the numerical simulation of dynamic soil-structure interaction problems both a frequency and a time domain formulation are presented. In order to be capable to consider more sophisticated models of the structure, the frequency domain algorithms for homogeneous and layered halfspaces have been coupled to the Finite Element Program ANSYS. Flexibility functions are presented for a concrete slab track system. Furthermore the stress distribution in the subsoil is calculated and visualized. The time domain formulation is applied for demonstrating the basic phenomena of a moving load passing by with sub- and supercritical speed. Besides that, a nonlinear, tension-free condition of contact between the track and the subsoil is mentioned briefly

Rapid GRB Follow-up with the 2-m Robotic Liverpool Telescope

We present the capabilities of the 2-m robotic Liverpool Telescope (LT), owned and operated by Liverpool John Moores University and situated at ORM, La Palma. Robotic control and scheduling of the LT make it especially powerful for observations in time domain astrophysics including: (i) rapid response to Targets of Opportunity: Gamma Ray Bursts, novae, supernovae, comets; (ii) monitoring of variable objects on timescales from seconds to years, and (iii) observations simultaneous or coordinated with other facilities, both ground-based and from space. Following a GRB alert from the Gamma Ray Observatories HETE-2, INTEGRAL and Swift we implement a special over-ride mode which enables observations to commence in about a minute after the alert, including optical and near infrared imaging and spectroscopy. In particular, the combination of aperture, site, instrumentation and rapid response (aided by its rapid slew and fully-opening enclosure) makes the LT excellently suited to help solving the mystery of the origin of optically dark GRBs, for the investigation of short bursts (which currently do not have any confirmed optical counterparts) and for early optical spectroscopy of the GRB phenomenon in general. We briefly describe the LT's key position in the RoboNet-1.0 network of robotic telescopes.Comment: 6 pages, 2 figures, to appear in the proceedings of Interacting Binaries: Accretion, Evolution and Outcomes, 4-10 July 2004, Cefalu, Sicily, Italy, eds. Antonelli et a

GMRT Observations of the 2006 outburst of the Nova RS Ophiuchi: First detection of emission at radio frequencies < 1.4 GHz

The first low radio frequency (<1.4 GHz) detection of the outburst of the recurrent nova RS Ophiuchi is presented in this letter. Radio emission was detected at 0.61 GHz on day 20 with a flux density of ~48 mJy and at 0.325 GHz on day 38 with a flux density of ~ 44 mJy. This is in contrast with the 1985 outburst when it was not detected at 0.327 GHz even on day 66. The emission at low radio frequencies is clearly non-thermal and is well-explained by a synchrotron spectrum of index alpha ~ -0.8 (S propto nu^alpha) suffering foreground absorption due to the pre-existing, ionized, warm, clumpy red giant wind. The absence of low frequency radio emission in 1985 and the earlier turn-on of the radio flux in the current outburst are interpreted as being due to higher foreground absorption in 1985 compared to that in 2006, suggesting that the overlying wind densities in 2006 are only ~30% of those in 1985.Comment: 14 pages, 1 figure. Accepted for publication in ApJ

Statistical kinetic treatment of relativistic binary collisions

In particle-based algorithms, the effect of binary collisions is commonly described in a statistical way, using Monte Carlo techniques. It is shown that, in the relativistic regime, stringent constraints should be considered on the sampling of particle pairs for collision, which are critical to ensure physically meaningful results, and that nonrelativistic sampling criteria (e.g., uniform random pairing) yield qualitatively wrong results, including equilibrium distributions that differ from the theoretical J\"uttner distribution. A general procedure for relativistically consistent algorithms is provided, and verified with three-dimensional Monte Carlo simulations, thus opening the way to the numerical exploration of the statistical properties of collisional relativistic systems.Comment: Accepted for publication as a Rapid Communication in Phys. Rev.

Hubble Space Telescope Imaging of the Expanding Nebular Remnant of the Recurrent Nova RS Ophiuchi (2006)

We report Hubble Space Telescope imaging obtained 155 days after the 2006 outburst of RS Ophiuchi. We detect extended emission in both [O III] and [Ne V] lines. In both lines, the remnant has a double ring structure. The E-W orientation and total extent of these structures (580+-50 AU at d=1.6kpc) is consistent with that expected due to expansion of emitting regions imaged earlier in the outburst at radio wavelengths. Expansion at high velocity appears to have been roughly constant in the E-W direction (v_{exp} = 3200+-300 km/s in the plane of the sky), with tentative evidence of deceleration N-S. We present a bipolar model of the remnant whose inclination is consistent with that of the central binary. The true expansion velocities of the polar components are then v = 5600+-1100 km/s. We suggest that the bipolar morphology of the remnant results from interaction of the outburst ejecta with a circumstellar medium that is significantly denser in the equatorial regions of the binary than at the poles. This is also consistent with observations of shock evolution in the X-ray and the possible presence of dust in the infrared. Furthermore, it is in line with models of the shaping of planetary nebulae with close binary central systems, and also with recent observations relating to the progenitors of Type Ia supernovae, for which recurrent novae are a proposed candidate. Our observations also reveal more extended structures to the S and E of the remnant whose possible origin is briefly discussed.Comment: 13 pages, 2 figures, accepted for publication in ApJ

Tree-Particle-Mesh: an adaptive, efficient, and parallel code for collisionless cosmological simulation

An improved implementation of an N-body code for simulating collisionless cosmological dynamics is presented. TPM (Tree-Particle-Mesh) combines the PM method on large scales with a tree code to handle particle-particle interactions at small separations. After the global PM forces are calculated, spatially distinct regions above a given density contrast are located; the tree code calculates the gravitational interactions inside these denser objects at higher spatial and temporal resolution. The new implementation includes individual particle time steps within trees, an improved treatment of tidal forces on trees, new criteria for higher force resolution and choice of time step, and parallel treatment of large trees. TPM is compared to P^3M and a tree code (GADGET) and is found to give equivalent results in significantly less time. The implementation is highly portable (requiring a Fortran compiler and MPI) and efficient on parallel machines. The source code can be found at http://astro.princeton.edu/~bode/TPM/Comment: 29 pages, includes 12 figures; to appear in ApJS. Source code at http://astro.princeton.edu/~bode/TPM

Double radio peak and non-thermal collimated ejecta in RS Ophiuchi following the 2006 outburst

We report Multi-Element Radio-Linked Interferometer Network, Very Large Array, One-Centimetre Radio Array, Very Long Baseline Array (VLBA), Effelsberg and Giant Metrewave Radio Telescope observations beginning 4.5 days after the discovery of RS Ophiuchi undergoing its 2006 recurrent nova outburst. Observations over the first 9 weeks are included, enabling us to follow spectral development throughout the three phases of the remnant development. We see dramatic brightening on days 4 to 7 at 6 GHz and an accompanying increase in other bands, particularly 1.46 GHz, consistent with transition from the initial ‘free expansion’ phase to the adiabatic expansion phase. This is complete by day 13 when the flux density at 5 GHz is apparently declining from an unexpectedly early maximum (compared with expectations from observations of the 1985 outburst). The flux density recovered to a second peak by approximately day 40, consistent with behaviour observed in 1985. At all times the spectral index is consistent with mixed non-thermal and thermal emission. The spectral indices are consistent with a non-thermal component at lower frequencies on all dates, and the spectral index changes show that the two components are clearly variable. The estimated extent of the emission at 22 GHz on day 59 is consistent with the extended east and west features seen at 1.7 GHz with the VLBA on day 63 being entirely non-thermal. We suggest a two-component model, consisting of a decelerating shell seen in mixed thermal and non-thermal emission plus faster bipolar ejecta generating the non-thermal emission, as seen in contemporaneous VLBA observations. Our estimated ejecta mass of 4 ± 2 × 10−7 M⊙ is consistent with a white dwarf (WD) mass of 1.4 M⊙. It may be that this ejecta mass estimate is a lower limit, in which case a lower WD mass would be consistent with the data